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Three-wave interactions form the basis of our understanding of many pattern forming systems 
because they encapsulate the most basic nonlinear interactions. In problems with two comparable 
length scales, it is possible for two waves of the shorter wavelength to interact with one wave of 
the longer, as well as for two waves of the longer wavelength to interact with one wave of the 
shorter. Consideration of both types of three-wave interactions can generically explain the presence 
of complex patterns and spatio-temporal chaos. Two length scales arise naturally in the Faraday 
wave experiment with multi-frequency forcing, and our results enable some previously unexplained 
experimental observations of spatio-temporal chaos to be interpreted in a new light. Our predictions 
are illustrated with numerical simulations of a model partial differential equation. 



PACS numbers: 47.54.-r, 47.52,+j, 05.45.-a, 47.35.Pq 

Patterns arise in many non-equilibrium physical, chem- 
ical and biological systems, often when a uniform state 
is subjected to external driving and becomes unstable to 
modes with a finite wavelength. In some systems, modes 
with a second wavelength can play an important role in 
pattern formation if these modes are either unstable or 
only weakly damped. The interaction between two waves 
of one wavelength with a third wave of the other wave- 
length is known both experimentally and theoretically to 
play a key role in producing a rich variety of interesting 
phenomena such as quasipatterns, superlattice patterns 
and spatio-temporal chaos (STC) [THE]- 

In this paper, we focus on three-wave interactions 
(3WIs) involving two comparable wavelengths and de- 
velop a criterion for when such interactions are likely to 
lead to STC, as opposed to steady patterns and quasi- 
patterns. The mechanism we describe is generic, and will 
apply to any system in which such 3WIs can occur, such 
as the Faraday wave experiment [JJf3], coupled Turing 
systems [I] and some optical systems [5]. In order to il- 
lustrate our criterion we focus on two examples. The first 
is the Faraday experiment, in which patterns of standing 
waves are excited on the surface of a fluid by periodi- 
cally forcing the fluids' container up and down. Using 
a multi-component forcing f(t) enables the excitation of 
waves with comparable wavelengths. Experimentally, the 
phases and amplitudes of the different components of f(t) 
have been shown to determine whether simple pattterns, 
superlattice patterns, quasipatterns or STC are seen [lj- 
[3]. A theoretical understanding of the stabilization of 
some superlattice patterns has been developed using a 
single 3WI [21 13 [7], but up until this point there has 
been no explanation for the presence of STC. We show 
how our extension of the notion of 3WIs can explain the 
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FIG. 1: Three- wave interactions between waves on two critical 
circles, with outer radius 1 and inner radius q > |. (a) A 
vector qi on the inner circle can be written as the sum of two 
vectors fci and h-z on the outer circle; this defines the angle 
6 Z = 2 arccos(g/2). (b) The vector fci is the sum of two inner 
vectors qi and qy, this defines the angle 9 W = 2 arccos(l/2g). 
(c) Similarly, &2 is the sum of two inner vectors, and each of 
these is in turn the sum of two outer vectors, and so on. 



origin of STC, quasipatterns and other features of the 
experimental results of [3]. Our second example is a new 
model PDE that has been designed to enable an explo- 
ration of 3WIs with two comparable wavelengths. 

We suppose that in the system two wavenumbers k — 1 
and k = q are excited with < q < 1. 3 Wis then can take 
two forms: first, two waves on the outer circle (wavenum- 
ber 1) interact nonlinearly with a wave on the inner circle 
(wavenumber q), for example fci and k 2 with q\ = ki+k 2 
(figure [lji) ; and second, if ^ < q < 1, two waves on the 
inner circle interact with one on the outer circle, for 
example vectors q 2 and q^ with k\ = q 2 + 53 (figure [I]d). 
The consequence of the presence of both forms of 
interaction when | < q < 1 is that, for almost all choices 
of q, nonlinear interactions necessarily generate infinitely 
many modes on both critical circles (figure [IJ2). This oc- 
curs because each vector on the outer circle is the sum of 
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two inner vectors, and in turn each of these inner vectors 
is the sum of two (different) outer vectors, and so on. 

So, for \ < q < 1, each mode on each circle will be 
influenced by two types of 3WI: one type involving in- 
teraction with a pair of modes from the other circle, and 
the second involving interaction with one mode from each 
circle. (There is also the usual interaction between three 
waves at 120° on a single circle.) To understand the effect 
of this multitude of nonlinear interactions (figure [lj:), we 
write the pattern-forming field U(x, y, t) as 

oo oo 

U= zA t Y kyX + J2 ^i(*)e i9ra5 + h.o.t., 

J=X,|fej 1=1 j = l,\ qj \=q 

where t is the time-scale for the slow evolution of the 
amplitudes, and the sums are over all modes on the two 
circles. In order to quantify the four coefficients that 
govern the two types of 3WIs, we write down archetypal 
equations (at lowest order) for the two different types of 
3WI, first involving a triad of wavevectors fci, k 2 and 
qi = k\ + k 2 , with amplitudes z±, Z2 and w\\ 

z\ = /izi + Q ZW Z2W\ + cubic terms, 
z 2 = jJiZi + Q Z wZ\W\ + cubic terms, (I) 
wi = vwi + Q zz z\Z2 + cubic terms, 

where [i and v are growth rates corresponding to wave- 
numbers I and q, and Q zw and Q zz are coefficients of 
the quadratic terms. Cubic terms play an important role 
in the dynamics of ([lj [5], and their effect is included 
in the discussion below. We write similar equations for 
the second type of 3WI, involving the triad q2, q% and 
fei = Q2 + Q3, with amplitudes W2, W3 and z\\ 

W2 = vw 2 + Q wz wj,zi + cubic terms, 
W3 = UW3 + Q wz w 2 zi + cubic terms, (2) 
z\ = fj,z% + Q WW W2W3 + cubic terms, 

where Q wz and Q ww are two more quadratic coefficients. 

These equations should be replicated for all possible 
combinations of modes from the two circles, whenever 
3WI is possible. The evolution of this infinite set of waves 
could be exceedingly complicated, but we argue it is likely 
to be most strongly influenced by the outcome of the two 
types of 3WIs governed by ([I]) and pj. To this end, we 
summarize the results of [S] and others on equations ([!]) 
with appropriate cubic terms. Pure mode solutions are 
stable for ranges of values of \i and v. The behavior of 
mixed solutions, where all three amplitudes are non-zero, 
is influenced heavily by the signs of the quadratic coeffi- 
cients. If these have the same sign (Q ZZ Q ZW > 0), there 
is a range of values of (/it, v) where steady solutions of ([!]) 
with all three modes non-zero are stable. In contrast, if 
the quadratic coefficients have opposite sign (Q ZZ Q ZW < 
0), there is a range of values of (/x, v) where there is time- 
dependent competition between the three amplitudes [9] . 

In a similar manner, in the multitude of 3WIs that oc- 
cur when I < q < 1, the behavior of each type of 3WI will 



depend on the products of the quadratic terms, Q ZZ Q ZW 
and Q ww Qw Z in fl} and (J2J, and will determine whether 
each type of interaction results in competition between 
modes or not. Consequently the signs of Q ZZ Q ZW and 
QwwQw Z will strongly influence the dynamics of the over- 
all system. Specifically, if Q ZZ Q ZW and Q ww Qw Z are both 
negative, then we expect to find time-dependent com- 
petition leading to intermittent appearance of patterns 
with 9 Z or 9 W , including the possibility of STC, with two 
full circles of modes present. Conversely, if Q ZZ Q ZW and 
QwwQw Z are both positive, only steady patterns are to be 
expected: these may be simple patterns such as stripes 
or hexagons, or they may have complex spatial struc- 
ture, with a superposition of waves having a variety of 
orientations featuring the special angles 9 W and 9 Z , pos- 
sibly with patches of one pattern in a background of the 
other. In the intermediate case where, Q ZZ Q ZW > and 
QwwQw Z < (or the other way around), then resonant 
triad interactions involving 9 Z are reinforced, while triads 
involving 9 W can be competitive (or vice- versa). In this 
case we expect to find steady superlattice patterns in- 
volving 9 Z , or time-dependent competition between such 
superlattice patterns with different orientations. 

We emphasise that these expectations for the overall 
behavior of the pattern-forming problem are based on the 
idea that the dynamics of any particular vector drawn 
from the infinite number of vectors on the two critical 
circles will be dominated by its linear growth rate and 
the 3WIs associated with the two types of triad to which 
that vector belongs. 

There are two special values of the radius ratio q that 
do not immediately generate an infinite number of waves: 
when q is g(\/5- 1) = 0.6180 or |(V6- V%) = 0.5176, 
corresponding to 9 Z = 29 w = 144° and 9 Z = 56 w = 150°, 
only a finite number of wavevectors on the critical circles 
is generated. These two special cases correspond to 10- 
fold and 12-fold quasipatterns respectively (although an 
infinite set of vectors is generated in the 10-fold case once 
hexagonal interactions are included). The values of q 
that would lead to 8-fold patterns generate an infinite 
number of vectors. Hexagons and squares are also not 
exceptions: for example, if q = l/y/2, we have 9 W = 90° 
but 9 Z = 138.6°, which does not fit in a square lattice. 

Considering both types of 3WIs provides an explana- 
tion of some of the observed Faraday wave phenomena in 
the experiments of [3]. These use the three- frequency 
forcing f(t) — a m cos (muit) + a n cos (mot + <f> n ) + 
a p cos (piot + (f> p ) , where (to, n,p) — (4, 5, 2) and (6, 7, 2). 
The 3WIs are between harmonic modes, driven by the 
(4, 2) and (6, 2) components of the forcing, with the pri- 
mary instability being to the larger frequency. The radius 
ratios are q ~ 0.52 and q 0.38 in the two cases: these 
are the radius ratios that arise in 12-fold quasipatterns 
and 22° superlattice patterns. The presence of these two 
patterns as a function of the phases (65, 4>2) and ((/>?, <fo) 
is indicated by light areas in figure [2^a,b). The slopes 
of the bands and their periodicity are related to time 
translation symmetries of components of f(t) [B]. 
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FIG. 2: (a,b) Experimental results of Ding and Umbanhowar [3], reproduced with permission, (a) The 30° angular autocorre- 
lation function with (4, 5, 2) forcing: white indicates 12-fold quasipatterns, as a function of (fa, fa). Black indicates disordered 
patterns, (b) The 22° angular autocorrelation function with (6, 7, 2) forcing: white indicates 22° superlattice patterns, as 
a function of (fa, fa). (c,d,e) Products of quadratic coefficients. (c,d) Parameter values as in the experiment in (a), with 
(c) QzzQzw (contours in steps of 0.2) and (d) Q WW Q WZ (contours in steps of 0.05) (e) Parameter values as in (b), showing 
QzzQzw (contours in steps of 0.1). In all cases, the thick line is the zero contour, solid (dashed) lines are positive (negative) 
contours. In (c,e), QzzQzw is negative only in very narrow diagonal bands; in (d), QwwQwz is positive in two horizontal bands. 



The two cases, (4, 5, 2) and (6, 7, 2) excitation, differ in 
that in the first case, q > \, and 'two-way' 3WIs between 
the different critical modes is possible, whereas in the sec- 
ond case, q < \, and no such two-way 3WIs can occur. 
We have calculated the values of the relevant quadratic 
coefficients using weakly nonlinear analysis of the Navier- 
Stokes equations for the experimental parameters by ex- 
tending the method in [TO]. In hgure [2]jc,d), we show 
QzzQzw and QwwQwz f° r the (4,5,2) case. Both prod- 
ucts can take either positive or negative values, although 
each is dominated by one sign. In figure [2](e), we show 
QzzQzw for the (6, 7, 2) case. Comparing with the exper- 
imental results in figure [2^a,b), we note that regions of 
12-fold quasipatterns and 22° superlattice patterns cor- 
relate extremely well with regions of positive QzzQzw- 
Conversely, when Q zz Qzw is small or negative, the auto- 
correlation function in the experiments is low and STC 
is seen [llj . The extra structure in figure [2ja), where the 
dark stripes broaden horizontally just above (\>i — ^ and 
np, is aligned with changes in QwwQwz in figure 2fd), 
although our expectation would be for enhanced sTC 
where QwwQwz is negative. 

There is clearly scope for further investigation of the 
role of the two-way 3WIs in this and other experiments, 
motivated by the encouraging alignment between the ex- 
perimental autocorrelation functions and the calculated 
quadratic coefficients. 

We have further investigated the importance of the 
signs of the quadratic coefficients in the two-way 3WIs 



by devising a new model PDE for a field U(x,y,t): 
— = C(p, v)U+Q 1 U 2 +Q 2 UV 2 U+Q 3 \WU\ 2 -U 3 . (3) 

The model ^ is an extension of those of [12] , modified 
to allow the growth rates of two critical modes to be 
controlled independently. The linear part of the PDE C 
acts on a mode e lkx with eigenvalue a(k). The depen- 
dence of the eigenvalue on k is specified by er(l) = \i and 
tj(q) = v, controlling the growth rates of the modes of 
interest; er'(l) = cr'(q) = 0; and er(0) = a < 0, control- 
ling the depth of the minimum between k = 1 and k = q. 
With a an even function of k, these requirements lead to 
a fourth-order polynomial in k 2 : 

<k) = k 2 (m^Bm + m j _ fc2)2) 

where A(k) = [k 2 (q 2 - 3) - 2q 2 + 4)(q 2 - k 2 ) 2 q 4 and 
B{k) = (k 2 {3q 2 - 1) + 2q 2 - 4q 4 )(l - k 2 ) 2 . The linear 
operator L is obtained by replacing k 2 by —V 2 . The non- 
linear terms in the PDE model are simple quadratic and 
cubic combinations of U and its derivatives. Standard 
weakly nonlinear theory gives the values of Qzz-, Qzw, 
Qww an d Qwz'- having the three quadratic terms in ^ 
enables different sign combinations to be chosen in the 
amplitude equations ([I]) and 

Simulations of equation pL carried out on a 30 x 30 
domain for a range of (//, v), are shown in figure [3]ja) 
for q — 0.66 and values of Qi, Qi and Qz that give 
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FIG. 3: (a) Bifurcation set showing the patterns that are 
seen for q — 0.66, (To = — 2, Q\ — 0.3, Qi = 1.3, Q3 = 1.7 (so 
QzzQzm < and QwwQwz < 0). We find z-hexagons (fc = 1), 
w-hexagons (k — q), spatio-temporal chaos (STC), mixed pat- 
terns (w-stripes with patches of z-rectangles) and two types 
of superlattice pattern: one with 6 modes on the outer circle 
and 12 on the inner, separated by W , and the other with 6 
modes on the inner circle and 12 on the outer, separated by 8 Z ■ 
Steady (time-dependent) patterns are indicated with a + (o). 
The (fi, v) scale is not uniform, (b) Pattern for parameters 
in the STC region, with fj, = v = 0.00707. The correlation 
length of this pattern is about 1-2 wavelengths, (c) The two 
critical circles are clearly seen in its power spectrum. 



QzzQzw < and QwwQwz < 0. We find a variety of 
patterns and STC, with a typical solution shown in fig- 
ure [3] (b,c). For other choices of q and the quadratic 



coefficients, the behavior of the PDE is broadly in line 
with the expectations. When q < |, we find no examples 
of STC. For q > |, we find steady and time-dependent 
patterns with many modes on both critical circles. In 
general, with both Q ZZ Q ZW and QwwQwz positive, these 
patterns are all steady; time dependence (chaos or STC) 
is much more common when both QzzQzw and QwwQwz 
are negative. We find steady and time-dependent quasi- 
patterns only with the special values q = 0.5176 and 
q = 0.6180. Further details will be presented elsewhere. 

We conclude that whenever there are 3WIs between 
waves on two critical circles with radius ratio between 1 
and 2, interactions in both directions must be taken into 
account. This is generic: these two-way 3 Wis will be 
important in any problem with two comparable length 
scales. Most values of q > h lead to the possibility of 
generating an infinite number of interacting waves. The 
exceptions are those values associated with 10- and 12- 
fold quasipatterns. The outcome of the 3WIs will be in- 
fluenced by the signs of the quadratic coefficients in Q 
and ([2]), with time-dependence and STC most likely in 
the case of both pairs of quadratic coefficients having op- 
posite sign. These ideas have been confirmed in numer- 
ical investigations of a model PDE pj). By computing 
quadratic coefficients from the Navier-Stokes equations, 
we have shown that these ideas are in broad agreement 
with the experimental results of [3]: as the phases in 
the forcing are varied, switches between 12-fold quasipat- 
terns or 22° superlattice patterns and STC line up with 
changes of sign of the products of quadratic coefficients. 
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